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A finite-length topological insulator nanowire, proximity-coupled to an ordinary bulk s-wave 
superconductor and subject to a longitudinal applied magnetic field, is shown to realize a one- 
dimensional topological superconductor with unpaired Majorana fermions localized at both ends. 
This situation occurs under a wide range of conditions and constitutes what is possibly the most 
easily accessible physical realization of the elusive Majorana particle in a solid-state system. 



Predicted in a seminal 1937 papei^^ as purely real (as 
opposed to complex- valued) solutions of the Dirac equa- 
tion describing a spin-^ particle, Majorana fermions are 
distinguished by the curious fact that they are their own 
anti-particles. More precisely, in second quantized for- 
mulation, the creation and annihilation operators for 
Majorana fermions coincide. In high-energy physics 
compelling theoretical arguments suggest that neutrinos 
might be Majorana fermions but a convincing experi- 
mental proof is yet to be given^l In condensed matter 
physics Majorana fermions can appear as emergent de- 
grees of freedom in certain systems of electrons when su- 
perconducting order or strong correlations are presentPEl. 
Over the past decade solid state realizations of Majorana 
fermions have been under intense theoretical study both 
as a fundamental intellectual challenge and as a possi- 
ble platform for fault-tolerant quantum computatiorP^. 
Yet their remarkable properties - and indeed their very 
existence - await experimental confirmation. 

The purpose of this Communication is to advance a 
proposal for a new type of solid-state device that can 
serve as a host for Majorana fermions under a wide range 
of experimentally accessible conditions. Our proposed 
device, depicted schematically in Fig. 1, draws con- 
ceptually on recent ideas to realize Majoranas in both 
two- and one-dimensional heterostructures composed of 
a topological insulator or a semiconductor with strong 
spin-orbit coupling^ coupled to an ordinary super- 
conductor through a proximity effect. Specifically, we 
consider a nanowire (i.e. quantum wire with nanometer- 
scale cross-section) fashioned out of a strong topological 
insulator (STI), such as Bi2Se3 or Bi2Te2Se, placed on 
top of an ordinary s-wave superconductor (SC), subject 
to applied magnetic field along the axis of the nanowire. 
Through a combination of analytical insights and numer- 
ical calculations we demonstrate below that when the 
magnetic fiux through the nanowire cross-section is close 
to a half-integer multiple of the fundamental flux quan- 
tum $0 = hc/e, the topologically protected surface state 
realizes a one-dimensional topological superconducto^^ 
with Majorana fermions localized near the ends of the 
wire. We note that single-crystalline Bi2Se3 nanowires 
with ribbon geometry (i.e. 'nanoribbons') have been 
fabricatecf^ and transport experiments in these show un- 
ambiguous evidence for the topologically protected sur- 
face states up to, possibly, a room temperatur^^^. Very 
recently, the SC proximity effect has been demonstrated 



in Sn-Bi2Se3 interface^^^. 

The device depicted in Fig. 1 appears superficially sim- 
ilar to t he de vices based on semiconductor wires proposed 
in Refs.liSllil. The physics underlying the emergence of 
Majorana fermions is nevertheless fundamentally differ- 
ent. In the semiconductor wires the applied magnetic 
field serves to open up a gap through the Zeemann cou- 
pling to electron spins whereas our proposal relies exclu- 
sively on the orbital effect of the applied field. This dif- 
ference leads to several key advantages of our device over 
the previous proposals. First, for the topological phase 
to occur in a semiconductor wire it is essential that the 
chemical potential is fine-tuned to lie inside the Zeemann 
gap, whose typical size in IT field is ~ 1 meV or less. In 
our device, by contrast, chemical potential can be any- 
where inside the TI bulk gap, which is ~ 300 meV in 
Bi2Se3. Second, as explained below, our device can be 
operated in the regime where the TI surface state is pro- 
tected by time reversal symmetry (T) and the induced 
pairing gap is therefore robust against non-magnetic dis- 
order. Such a pr otection is absent in ordinary semicon- 
ductor wiresISIIi]. Finally, the energy gap protecting the 
Majorana fermions in our setup is generically an order of 
magnitude larger than in previous proposals. 

We begin with a qualitative discussion of the physics 
behi nd th e proposed device. As explained in previous 
workJ^H the key ingredient necessary to build a topolog- 
ical superconductor is an underlying normal state char- 
acterized by the electron dispersion with the spin degen- 
eracy removed. It is possible to achieve this situation 
on the surface of a TpSHIS. The low-energy fermionic 
excitations on such a surface are governed by the Dirac 
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FIG. 1: Schematic of the proposed device. Magnetic field B 
is applied along the axis of the wire taken to coincide with 
the z-direction. 
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FIG. 2: a) Surface state excitation spectra Eti for various val- 
ues of magnetic flux $ — 7?$o- Solid and dashed lines indicate 
doubly degenerate and non-degenerate bands, respectively, b) 
Kitaev's Majorana number: A4 = —1(4-1) in shaded (white) 
regions. Numerals inside the squares indicate the number of 
Fermi points for fc > 0. c) A possible shape of SC/magnetic 
domain wall located near z = 0. Dashed line shows the exact 
zero-mode solution u{z) for this domain wall. 



Hamiltoniarf^ 

h = ^w[/iV • n + n • (p X s) + (p X s) • n] -|-s • m, (1) 

where p = — ifiV, v is the Dirac velocity, n unit vector 
normal to the surface, and s — (si,S2,S3) the vector of 
Pauli spin matrices, m denotes the magnetization vec- 
tor, caused e.g. by the Zeemann coupling of spins to the 
external magnetic field. Inclusion of the latter is not es- 
sential for the functionality of the proposed device but 
will prove useful in subsequent considerations. 

Now consider a TI wire in the shape of a cylinder with 
radius R and magnetic field B applied along its axis. The 
magnetic flux is included by replacing the momentum 
operator with tt — p — (e/c)A, where A = 77$o(^ x 
r)/27rr^ is the vector potential. $ = 77<I>o represents the 
total magnetic flux through the cylinder. Taking fi = 
(cos (/3, sin t/?, 0) and m = the spectrum of Hamiltonian 
Q readPl 



Eu = ±vh\ 



i?2 



(2) 



Here k labels momentum eigenstates along the cylinder 
while I = 0, ±1, . . . is the angular momentum. The spec- 
trum in Eq. ([2| is clearly periodic in 77 which reflects the 
expected $o-periodicity in the total flux. Our identifi- 
cation of the suitable 'spinless' normal state hinges on 
the following observation. For = all branches of E)^i 
are doubly degenerate (Fig. 2a). For 7^ 0, however. 



the degeneracy is lifted and one can always find a value 
of the chemical potential fi that yields a single pair of 
non-degenerate Fermi points, as illustrated in Fig. 2a. 
Pairing induced by the proximity effect in such a state 
is then expected to drive the system into a topological 
phase. 

One can formalize the above argument by considering 
Kitaev's Majorana number 7W defined as = (— 1)", 
where v represents the number of Fermi points for fc > 0. 
In the limit of weak pairing, = — 1 indicates the exis- 
tence of unpaired Majorana fermions at the ends of the 
wir^i^. Fig. 2b shows M. calculated from the spectrum 
Eq. ^ as a function of ^ and 77. We observe, specifically, 
that when 77 = 1/2, i.e for the flux equal to half-integer 
multiple of $0, Majorana fermions will appear for any 
value of the chemical potential (as long as it lies inside 
the bulk gap). This result is easily understood by not- 
ing that for ?7 = 1/2 the gapless I = branch is non- 
degenerate while the remaining branches are all doubly 
degenerate. Thus, the number of Fermi points for fc > 
is odd for any value of /x. It is also worth noting that, for 
the surface state, 77 = 1/2 represents a T- invariant point 
and the above pattern of degeneracies should therefore 
be robust with respect to non-magnetic disordeil^i^f^Sl jjj 
this situation Cooper pairs are formed from time-reversed 
electron states and the pairing gap is protected against 
disorder by Anderson's theorem. Below, we will explicitly 
demonstrate the existence and the robustness of the Ma- 
jorana fermions both analytically within the low-energy 
theory based on Hamiltonian ([T]) and numerically using 
a minimal lattice model. 

Writing the Hamiltonian ([I]) in cylindrical coordinates 
and with the ansatz for the wavefunction 
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the spinor ipki — {fki,9ki)'^ is an eigenstate of 



S2k + S3[{1 



rj) /R + m] 



(3) 



(4) 



Here we take v = h = 1 and m = mz. To illustrate the 
emergence of Majorana fermions in the simplest possible 
setting we now focus on the 77 = 1/2 case and consider 
chemical potential |^| < vh/R, i.e. intersecting only the 
I ~ branch of the spectrum Eq. (pi). The Hamiltonian 
for this branch then becomes hk = (ks2— fJ,)+rns3, where 
we have explicitly included the chemical potential term. 

With this preparation we can now construct 
the Bogoliubov-de Gennes Hamiltonian describing 
the proximity-induced superconducting order in the 
nanowire. In the second-quantized notation it reads 



H 



Ek H-^k^k 



ith = ifk^gkjlk^gikV and 



-Hk = 



hk 
-A* 



-h*, 



(5) 



In the following we consider the simplest s-wave pair po- 
tential Afc = Aois2 with Aq a (complex) constant or- 
der parameter, which corresponds to the pairing term 
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^o{flg-k~ 9kf-k) ■ Introducing Pauli matrices Tq, in the 
Nambu space we can write, assuming Aq real, 

"Hfc = r3(s2fc - + - r2S2Ao, (6) 

^2 + A2±2(fcV + 



with eigenvalues Ek — ±[k'^ ■ 
Ai2m2 + m2A2)i/2]i/2_ 

In the special case when fi — the spectrum simplifies. 



(7) 



The form of the spectrum above suggests that a local- 
ized zero-mode will exist at a boundary between SC and 
magnetic domains, i.e. when (m ± Aq) changes sign. 
We thus seek a zero-energy eigenstate H'ifo{z) = of 
v. = T3[-S2idz + S3m{z)] - r2S2A(z) with m(z), A(z) of 
the form indicated in Fig. 2c. A single Jackiw-Rossi 
zero mode'^' indeed exists and has the form '^o{z) ~ 
(1, — 1, 1, — l)-^u(z) with u{z) = uoexp dz'[A{z') ~ 
m(z')] and uq a normalization constant. The associated 
field operator 



u{z)[f{z)-giz) + f\z)-gHz)]dz, (8) 



has the property ipQ — 4'o a-nd represents, therefore, a 
Majorana fermion. 

The above explicit calculation establishes the existence 
of an unpaired Majorana mode at a SC/magnetic domain 
wall in a TI nanowire under very special conditions. We 
now argue that the effect is in fact generic. First, we rea- 
son that the magnetic order, although convenient in the 
derivation, is in fact irrelevant. Consider a nanowire of 
length L ^ ^, the lengthscale of the zero mode, with the 
domain wall located near its center. In a physical system 
Majoranas always come in pairs. Since the second Ma- 
jorana fermion evidently cannot live in the gapped bulk 
(or at the magnetic end) we conclude that it must be 
localized at the SC end, irrespective of the exact bound- 
ary condition. Second, it is easy to see that the chemical 
potential can be moved away from zero without perturb- 
ing the Majoranas. Indeed, with m = the spectrum of 
Eq. ^ reads Ek = ±[(fe ± fi)^ + AgJ^/^ indicating that 
the bulk of the wire shows a SC gap Aq for any value 
of /i. Therefore, Majorana end-states will persist even as 
fi is varied away from 0. When the chemical potential 
intersects additional bands then each band contributes 
a single Majorana end-state. Any even number of these 
will pair up to form ordinary fermions (whose energies 
will generically be non-zero) but for an odd number of 
occupied bands a single unpaired zero-energy Majorana 
will remain. This consideration elucidates the physical 
meaning of Kitaev's Majorana number A4. We note that 
the topological nature of the zero mode guarantees its 
stability against smooth deformations of the nanowire 
shape, as long as its bulk remains gapped and the total 
magnetic flux seen by the surface state is unchanged. 

To explicitly address the existence and robustness of 
Majorana end-states we now study the nanowire using 
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FIG. 3: a) Energy dispersion for an infinitely long TI wire 
with 20 X 20 base described by lattice Hamiltonian ([9| in the 
normal state with 77 = 0.52. For clarity only the low-energy 
portion of the spectrum is displayed in a part of the Bril- 
louin zone. Inset shows the spin expectation values for the 
gapless state at small positive k. The length of the arrow is 
proportional to the wavefunction amplitude, b) Lines sepa- 
rating regions with different Majorana number A4 — — 1(+1) 
extracted from the spectrum. All energies are in units of 
A = 150meV and we use parameters = t = 1, e = 4 and 
g — 32, corresponding to the strong TI phase with Z2 index 
(1;000) and bulk bandgap 2 A = 300meV. 



a concrete lattice model of Bi2Se3 family of materialJ^^. 
Specifically, we use the model given by Fu and Ber^^Il 
regularized on a simple cubic lattice, defined by a /c-space 
Hamiltonian 

= A//k(Ti-|-A(T3(s2sinA:j;-sisinfcj^)-|-A2CT2sinfc^, (9) 

with Mk = £ — '2tJ2a cos ka- Here ctq, represent the Pauli 
matrices acting in the space of two independent orbitals 
per lattice site. For A, > and 2t < e < 6t the 
system described by Hamiltonian ([o]) is a TI in Z2 class 
(I;000), i.e. a strong topological insulator. The magnetic 
field enters through the Peierls substitution, replacing all 
hopping amplitudes as tij — >■ tij exp [— (27rj/<I>o) // A • dl] 
and the Zeemann term — g/iBB-s/2 where fiB — eh/2mf.c 
is the Bohr magneton. In the SC state the BdG Hamil- 
tonian takes the form of Eq. ^ with At = Aois2 de- 
scribing on-site spin singlet pairing. 

We have solved the problem posed by Hamiltonian ^ 
in various wire geometries by exact numerical diagonal- 
ization and by sparse matrix techniques. Fig. 3a shows 
a typical example of the excitation spectrum in an in- 
finitely long wire with WxW cross-section in the normal 
state. We observe that for rj close to 1/2 the surface state 
is indeed gapless and the low-energy modes exhibit the 
expected pattern of degeneracy. Because of the surface 
state penetration into the TI bulk the surface electrons 
see a slightly smaller magnetic flux than the nominal flux 
$ = BW^ given by the wire geometry and the gapless 
state is shifted to a slightly higher value of rj. This is 
also seen in Fig. 3b which displays the Majorana num- 
ber for the same system. This figure indicates that for 
r] = 0.52(2n -I- 1) the system will be a ID topological SC 
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FIG. 4: a) Energy bands for an infinitely long TI wire with 
6x6 base in the SC state with rj = 0.49, fj. = 0.09, Ao = 0.08 
and g = (solid lines), and the energy levels for L = 36 
finite-length wire with open boundary conditions (red circles) 
obtained by exact numerical diagonalization. b) Three lowest 
positive energy eigenvalues obtained by the Lanczos method 
as a function of L. Dashed line represents the envelope func- 
tion O.OSQe"-^/^ with ^ = 17.5. 

for any value of fi inside the bulk gap. 

Superconducting order opens a gap in the electron ex- 
citation spectrum as illustrated in Fig. 4a. For an open- 
ended wire, crucially, our calculations reveal a pair of 
non-degenerate states at ±Eo inside the SC gap whose 
energies approach zero for large L as i?o oc e~^l^ . Fig. 
4b illustrates this exponential decay (which is in addi- 
tion modulated by oscillations at 2kp\ Higher energy 
eigenstates approach non-zero values close to Aq for large 
L. We have verified that the appropriate linear super- 
positions of the wavefunctions associated with the ±i?o 
eigenvalues are exponentially localized near the ends of 
the wire. The corresponding field operators then satisfy 
the Majorana condition i/i^^ — tp, and represent, up to 
exponentially small corrections in their separation, the 
Majorana zero modes. 

We conclude with comments on the experi ment al re- 
alization. For the existing Bi2Se3 nanowireJ^^'i^ with 
cross-sectional area 5 « 6 x 10~^^m^ the surface level 
spacmg IS SEs ~ vh^VfS ~ 7meV. At half flux quan- 



tum, which corresponds to the magnetic field strength 
B = $0/25 ^ 0.34T, the Zeemann energy scale iEz = 
girh? /2mf,S ~ 0.6meV (taking w 5 x lO'^m/s and 
g = 32) and thus probably negligible. Experiments 
on planar Sn-Bi2Se3 interfaces'^ show induced SC gap 
^ 0.2meV, a significant fraction of the native Sn bulk gap 
(^ 0.6meV). It thus appears conceivable that a pairing 
gap of several meV could be induced in Bi2Se3 nanowires 
by using in place of Sn a superconductor with larger 
bulk gap, such as NbTiN or MgB2. A gap of this size 
should permit detection of the Majorana fermion by scan- 
ning tunnelling techniques. The latter will show up as 
a zero-bias peak exponentially localized near the end of 
the wire in the topological phase, but will disappear as 
the phase boundary into non-topological phase is tra- 
versed by tuning either chemical potential or magnetic 
field. Unambiguous detection of Majorana fermions can 
be achieved by probing 47r-periodic Josephson current 
through a weak link in the wire as described in Refs.'*' 
Another intriguing possibility is to fabricate nanowires 
from Cua;Bi2Se3 which becomes a superconductor below 
4K[2S while simultaneously retaining the protected sur- 
face state^^ZI, 

Energy scales order of magnitude higher compared to 
the ordinary semiconductor wires, significantly reduced 
requirements for the chemical potential control and the 
sample purity (afforded by operation in the T-invariant 
regime) make TI nanowires promising candidates for the 
first experimental detection of Majorana fermions. Our 
numerical simulations of the model Hamiltonian ^ with 
disorder (to be reported separately) show essential ro- 
bustness against non-magnetic impuritie^^ll but system- 
atic studies of the effects of disorder and interactions 
along the lines of recent work a^^ ' ^" ' constitute an exciting 
future research direction. 
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